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Abstract 

Recent experiments have demonstrated that highly charged ions can be guided through insulating 
nanocapillaries along the direction of the capillary axis for a surprisingly wide range of injection 
angles. Even more surprisingly, the transmitted particles remain predominantly in their initial 
charge state, thus opening the pathway to the construction of novel ion-optical elements without 
electric feedthroughs. We present a theoretical treatment of this self-organized guiding process. 
We develop a classical trajectory transport theory that relates the microscopic charge- up with 
macroscopic material properties. Transmission coefficients, angular spread of transmitted particles, 
and discharge characteristics of the target are investigated. Partial agreement with experiment is 
found. 

PACS numbers: 34.50.Dy 
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I. INTRODUCTION 



The transmission of multiply and highly charged ions (HCI) through nanocapillaries has 
recently been employed as a tool to study the interaction with surfaces, specifically, the 
internal walls of the capillary. The original motivation was to delineate the initial state 
of hollow-atom formation at large distances from the surface information hardly 

accessible by scattering at fiat surfaces. Image attraction and close collisions tend to erase 
the memory on this early state of ion-surface interaction when the ion suffers close encounter 
as it is either reflected from the topmost layer of the surface or penetrates into the bulk. A 
previous first step to access information on the early stages of the neutralization process was 
the usage of "stepped", or terrace-decorated, surfaces 41]. Capillaries provided an attractive 
alternative as the ion can escape the capillary prior to a close encounter with the surface, 
the hollow atom can be 

focused on metallic capillaries with typical radii of r = 50 to 100 nm ("nanocapillaries") 
and a length L of about one /im, thus featuring an aspect ratio of the order of 1:20 and 
geometric opening angles 6q < 3°. Incidence angles 6in larger than 9q therefore necessarily 
lead to the impact of all projectiles on the surface and, consequently, to the destruction of 
the hollow atom (or hollow ion). 

More recently, capillaries through insulating foils (PET or "Mylar" 6|1) and Si02 [7] with 



aspect ratios ~ 1 : 100 have been studied in several laboratories ja |9|, |lO| . Unexpectedly, 
considerable transmission probabilities for projectiles in their initial charge state were mea- 
sured for incidence angles as large as ~ 20°. Apparently, ions are guided along the capillary 
axis with a spread (FWHM) of AOout of several degrees for mylar y] but close to geometric 



opening 9o for Si02 j8[. Keeping the initial charge state, contrary to the expected neutral- 
ization upon approach of the internal capillary surface, suggests that the ions bounce off the 
walls at distances larger than the critical distance Rc ~ y/2Q/W \'^: {Q: charge state; W: 
workfunction of capillary). Key to this process is the charging up of the internal insulator 
walls due to preceding ion impacts. Ion guiding through the capillary ensues as soon as a 
dynamical equilibrium of self-organized charge up by the ion beam, charge relaxation, and 
reflection is established. 

A theoretical description and simulation of this process poses a considerable challenge in 
view of the widely disparate time scales simultaneously present in this problem: 
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1. The microscopic charge- up and charge hopping due to the impact of individual ion 
impact takes place on a time scale of sub-/s to fs with typical hopping time Th < 10~^^ 
s. 

2. The transmission time Tt of a projectile ion through the capillary for typical ion energies 
of ^ 200 eV/u is of the order of n ^ 10^^° s. 

3. Typical average time intervals At between two subsequent transmission (or impact) 
events in the same capillary are, for present experimental current densities of nA/mm^ 
of the order of At ^ 0.1 s. 

4. Characteristic (bulk) discharge times for these highly insulating materials, can be 



estimated from conductivity data 
reach days. 



Ill ] to typically exceed > 10^ s and can even 



As this multi-scale problem spans a remarkable 18 orders of magnitude, a fully microscopic 
ab initio simulation covering all relevant scales is undoubtedly out of reach. The more modest 
goal of the present approach is therefore to develop a simulation that allows to interrelate 
the microscopic description of ion-surface impact with macroscopic properties of charge-up 
and transport. Aim is therefore to employ in our simulation only parameters that are not 
freely adjustable but can, at least as far as their order magnitude is concerned, be deduced 
from data for macroscopic material properties of the nanocapillary material. Specifically, 
the bulk discharge time as well surface charge diffusion constant will be estimated 
from surface and bulk conductivity data for mylar IjJ . 

The present approach represent a mean-field classical transport theory [12j based on a 
microscopic classical-trajectory Monte Carlo (CTMC) simulation for the ion transported, 
self-consistently coupled to the charge-up of and charge diffusion near the internal capillary 
walls. As a limiting case we also consider a simplified rate-equation approach. We find 
partial agreement with experimental data. We conclude this paper by pointing to future 
experimental and theoretical challenges to be overcome. 



3 



FIG. 1: lon-nanocapillary interaction, schematically: array of nano capillaries oriented along the 
surface normal with close-up of an individual capillary. An insulating material (PET) with di- 
electric constant e is covered on both sides with gold layers (dark shaded) preventing charge up 
of the target during experiment. Capillaries with radius a = 50 nm and L = 10 /im are typically 
D = 500 nm apart. Projectiles enter and exit the capillary under angles Om and 6 out with respect 
to the capillary axis. The capillary axis is either normal to the surface or Gaussian distributed 
with Ma < 2° (FWHM). 

II. METHOD OF SIMULATION 
A. The Scenario 

The simulations presented below incorporate the features of the following underlying 
scenario: an ensemble of ions of charge state Q are incident relative to the surface normal 
of the nanocapillary foil with angle 9 in. The axis of the capillary with radius a is assumed 
to be either aligned with the surface normal (Fig. [1} or Gaussian distributed around the 
orientation of the normal with width A^^. The lateral inter-capillary spacing D is of the 
order of -D ~ 10 a such that, on an atomic level, inter-capillary interaction effects can be 
neglected. Mesoscopic field effects resulting from the charge-up of the ensemble of capillaries 
will, however, be taken into account. The entrance and exit surfaces are covered by a 30 
nm metallic gold layer in order to avoid macroscopic charge-up of the capillary foil. The 
insulating material on the inside is characterized by a dielectric constant e = 3.3. In view 
of the low projectile velocities Vp for which guiding is studied, fp < 0.1 a.u., the static limit 
e{uj ^ 0) = e of the dielectric function e{uj) is appropriate to characterize the linear response. 
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About 4 % of the entrance surface is covered by capillary openings. Typical current densities 
are ~ nA/mm^ and the beam spot size is of the order of ~ mm^ simultaneously illuminating 
~ 10^ adjacent nanocapillaries. Ions entering the nanocapillary will be attracted by the 
image field 

^"^ A\(P\e + l ' ^' 
where d is the distance to the capillary wall. Therefore, even ions entering the capillary with 

angles 9in smaller than the geometric aspect ratio A, 

tan ^ 6*0 = A = a/L (2) 

will eventually experience close encounters with the surface. Therefore, the effective angle for 
transmission would be smaller than given by Eq.|21 unless guiding effects become operational. 
Ions reaching the surface will undergo a charge transfer sequence yielding a neutralized 
projectile and deposite a number Qeff of positive charges with Qeff ^ <5 at the surface. 
In our simulation, we use Qeff = Q for the microscopic charge up per ion impact at the 
capillary wall above all because the secondary electron emission coefficient 7 for impact of 

n 

slow ions at insulator surfaces is small 7^1 ^9']. We have, however, tested the influence 
of larger 7 by performing simulations in which up to 2Q positive charges and Q electrons 
with energies up to 50 eV were produced. Trajectories of the latter were followed and these 
electrons were allowed to neutralize positive charges on the capillary wall. No significant 
differences between simulations for different 7 could be found. 

The charged-up micro-patch will undergo surface diffusion with diffusion constant Dg as 
well diffusion into the bulk with diffusion constant Db. Bulk diffusion is extremely slow 
for highly insulating materials while the surface diffusion towards the grounded metallic 
layers will be a factor ~ 100 faster, thus governing the overall discharge process. Self- 
organized guiding sets in when a dynamical equilibrium between charge-up by a series of 
ion impacts at internal walls with an average frequency (At)^^ ^ 5 Hz and the charge 
diffusion is established such that the electrostatic repulsion prevents further impacts. The 
ion is reflected at distances from the wall larger than the critical distances from the surface 

Rc = V^/W . (3) 

The wall forms then an effective mesoscopic "trampoline" for subsequent ions and guides 
the projectile towards the exit. 



B. Projectile Trajectories 



We consider the propagation of 1 < n < ions through 1 < m < M nanocapillaries. 
The calculation of the n*^ ionic trajectory through the m*^ capillary proceeds by solving 
Hamilton's equation of motion 

f)J-f 

i^^) = (5) 
with random initial conditions subject to the energy constraint 

The non-trivial aspect of this classical-trajectory Monte Carlo simulation is that memory 
effects must be built into the potential V (Eq. The force experienced by the n'^^ ion 
depends on the history of all previous trajectories {R}^^} through the same (m = m') 
capillary as well as on the ensemble of neighboring charged-up capillaries m ^ m' via 
the mean field Fmean- It furthermore depends on the surface and bulk diffusion constants 
Dg^b of the material. By deducing values for Dg^b from macroscopic conductivity data and 
considering them to be independent from the stage of charge-up, we treat the discharge 
process in linear response, i.e. we neglect non-linear processes discussed recently 14, Ifij ]. 
Details of the determination of V will be discussed below. 

Solving Eqs. |3] and El is, while straight-forward in principle, quite demanding in view 
of the disparate time scales and ensemble sizes involved. We have therefore developed an 
algorithm that approaches the accuracy of conventional Runge-Kutta methods while con- 
siderably increasing the computing efficiency. First, the capillary is divided into three-sided 
prisms which are further subdivided into three tetrahedrons (see Fig. 12)). The force (Eq. Ej) is 
calculated before the simulation of a given trajectory only at the corners of the tetrahedrons. 
The starting point of the trajectory is determined from a uniform probability distribution 
over the entrance surface of the capillary, the entrance angle 6'j„ is selected from a Gaussian 
distribution centered around the nominal entrance angle. Its full width at half maximum 
(FWHM A,„ . 0.5", wa. U.en fro. expenn,enta, fl. TKe p.opa,a«on . 

calculated in three steps: first, the point of entrance into a tetrahedron is determined. Next, 
the average electric field calculated from the potential (Eq. in this structure is deter- 
mined and, as the last step, the exit point is derived from the analytical propagation of 



U(x1) 
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FIG. 2: Algorithm for propagation through the capiUary, schematicahy. The capillary is divided 
into tetrahedrons through which the trajectory is propagated analytically (cf. text). 

the trajectory within the tetrahedron. This exit point becomes the starting point of the 
trajectory within the next tetrahedron and so forth (see Fig. |21). From comparison between 
our method and a fourth-order Runge-Kutta simulation we have determined the optimal 
size of the tetrahedrons. We find a significant speed-up (up to a factor 100) while preserving 
accuracy to a good degree of approximation. 

The propagation of the trajectory is stopped if one of the following two requirements is 
met: Either the projectile reaches the exit surface of the capillary at which point the exit 
angle for this trajectory is calculated and the transmission probability of the ensemble is 
updated or the trajectory is not sufficiently deflected by the field of surface charges and 
hits the surface. The point of impact is determined and serves as starting point of the 
simultaneously performed calculation of the charge diffusion. Diffusion of the deposited 
Q charges is determined for a time interval At between two the subsequent trajectories 
(n, n + 1) enter the same capillary m. We have performed simulations both for Poissonian 
distributed intervals with mean At as well as for fixed intervals At = At. The results are 
indistinguishable which is not surprising as At is 9 orders of magnitude larger than Tf. 

C. Calculations of the Fields Inside the Capillaries 

The forces governing Eq. |S1 are determined from the electric fields F resulting from the 
charge- up as well as the image field. 




(7) 
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with given by Eq. [Tl Fcoui represents the electric field due to the charge-up of the wall 
of the rn}^ capillary, and Fmean is the mesoscopic field resulting from the charge-up of the 
entire ensemble of capillaries. 

Calculation of the charge distribution makes use of the macroscopic properties of these 
highly .usulatmg materials Q. Startmg point is the bulk couduot.vity ofMylar ^ 
10 ^^Vt ^m ^. The surface conductivity as is larger by a factor of about 100 |ll|. From the 
Einstein relation 

2 

ne , 
(^b,s = -j^Db,s (8) 

and using the experimentally determined value for the amount of charge carriers in Mylar n = 
lQi8 jj^ should be of the order of 10"^'^ m^s"^ and Dg ^ 10"^^ m^s~^ We use Df, = 

2 ■ 10~^^ m^s~^ and Dg = 100 ■ Db unless otherwise stated. The disparity of surface and bulk 
diffusion will be key for the understanding of the charge-up and discharge characteristics. 
The point to be noted is that for such low conductivities and large uncertainties in the 
effective number of carriers present, Eq. |Hlcan only serve for an order of magnitude estimate. 
We will later check on the consistency of these estimates by comparison with results of our 
simulation. After deposition of Q charges per impact, charges are assumed to undergo an 
unbiased random walk resulting in the diffusive spread of the charges. For simplicity and to 
avoid the introduction of any free parameter, modification of the random walk by mutual 
Coulomb repulsion is neglected. Accordingly, the probability distribution for the charge 



distribution in the surface a time interval At after impact is given by 



where p is the distance on the surface. Likewise, the distribution for bulk diffusion in 3D is 
given by 

with u the distance in 3D. Since surface diffusion is much faster and effectively controls the 
discharging at the boundary to the metallic layer, we treat surface diffusion explicitly using 
Eq. ini By contrast, we simulate bulk diffusion after the random walk into the bulk by a 
switch-off of the charge with a time constant Tf, estimated from as 

- - 1 (") 
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where we use as characteristic hopping length into the bulk half of the mean inter-capillary 

n 

distance (i.e. li, = D/2 ^ 5a 6]). 

The wall potential seen by the projectile entering the capillary is expressed by 

V.Mr-,t)= [ rfa;g4+^ "^"P'-<'-.';'''"f (12) 

f ij-t ij-ti ^^^^^ ly 'y ■ \ 



surface 
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In the first term of Eq. the surface charge density is given in terms of the sum over all 
surface diffusion distributions resulting from impacts at times tk < t, 

a(f,t) = (l-P,) J2 Psir-fk,t-tk). (13) 

k 

{tk<t) 

The second term in Eq. describes the exponentially decaying Coulomb interaction of the 
fraction of those {]} charges disappearing into the bulk at time tj and at position Vj with 
probability Ph- From Eq. [T21 the field Fyjaii = —^Vwaii can be determined. 

Finally, the mean field Fmean takes into account collective field effects of an entire ensemble 
of nanocapillaries approximately ~ 10^ of which are located within the beam spot. As they 
get simultaneously charged-up, they generate a mesoscopic electric field that is aligned along 
the capillary axis. Due to the large distance of the charge patches near the entry side from 
the exit surface (~ 10 yU m) the ensemble can be viewed as a charged condenser consisting of 
two metallic plates (the gold layers) with a dielectric material characterized by the dielectric 
constant e in between (Fig. Q). The mean field can be estimated from the total charge 
deposited on the capillary walls, the dielectric constant of the insulating material, and the 
irradiated target area. 

— * 

We have performed simulations of the potential and field Fmean near the exit surface 
using the program POISSON/SUPERFISH [l^. Equipotential lines from these calculations 
are shown in Fig. El For efficient implementation into the classical transport simulation 
we also developed a simplified analytic model for Fmean- To this end, the exit plane of 
the capillary is approximated by a infinitely thin metal plate with a hole of radius a. For 
2; ^ — oo, a homogeneous field Fq proportional to the amount of charge deposited at the 
entrance surface along the capillary axis is imposed as boundary condition, near the exit 
surface (z = 0) Dirichlet boundary conditions are applied. Then, the potential in the half 
sphere with 2; > is given by |3] 

Fn 

$+(r, 2) = -\h- z- arctan(6/^)] (14) 



FIG. 3: A metal plate with a hole with radius a separates two half spaces. On the right hand 
side {z > 0) Dirichlet boundary conditions are assumed for r, z — > oo. The boundary condition for 
z — > — oo is given by the requirement that the electric field converges to a homogeneous field -Fo||i. 



with 



Vc 
1 



The potential in the negative half space is determined by 

$"(r, -z) = $+(r, z) - FqZ 



with the matching conditions ai z = 



dz 



9$H 



dz 



[r < a). 



(15) 
(16) 

(17) 

(18) 



The analytic model gives fields Fmean = — V$ in close agreement with the numerically 
calculated field distribution. We therefore employ Eqs. El^Elin our CTT. The importance 
of the inhomogeneity of Fmean lies in its influence on the angular distribution of the guided 
ions by defocusing the transmitted beam near the exit surface. 

III. RESULTS 



A. Transmission function for ion guiding 

Following the pioneering work of Stolterfoht et al. several groups studied the ion 
guiding through Mylar (PET) 0, Q and Si02 nano capillaries j?!. Key finding is the 
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a) 



b) 



FIG. 4: Scatter plot of deposited charges in the interior of an individual capillary and typical 
trajectories for Om = 3°. a) zig-zag distribution leading to blocking (for Ds = Db); b) patch 
distribution leading to transmission (for Ds = WODh). Note that in b) two groups of typical 
trajectories can be observed, leading to two maxima in the angular distribution (Fig. Ions 
enter from the left, aspect ratio a/L = 100, not to scale. 

build-up of a self-organized charge distribution that enables the collision-free transmission. 
Patterns of the self-organized charge distribution are shown in Fig. |3] While Fig. 0^ shows a 
typical charge-up condition leading to blocking (electrostatic "bottleneck"), Fig. displays 
a distribution giving rise to guiding. The bottleneck condition was realized by setting Dg 
to be equal to the bulk values of Db, corresponding to a slow overall discharge time set here 
to be Tb ^ 35 min. In fact, such a slow discharge time has been suggested by experiment 
0] for the time interval over which recharging after beam switch-off could be realized. Only 
in the simultaneous presence of both fast (via surface transport, D^) and slow decay (via 
bulk transport Db) can guiding be established and maintained. When taking into account 

)oth surface and bulk diffusion, assumptions about a non-linear discharging characteristics 

l^l 14 , ^ need not to be invoked. 

The resulting transmission probability as a function of the angle of incidence relative 
to the capillary axis (Fig. Ej) displays significant transmission for angles well outside the 
geometric opening angle ^o- The quantitative agreement with the experimental transmission 
function is reasonably good considering the discrepancies between different data sets. The 
latter is, in part, due to the extraordinarily long bulk discharge times (from Db estimated 
to be in the range of hours) which makes measurements under reproducible conditions of 
complete initial discharge a challenge. 
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FIG. 5: Transmission function (transmission probability) as a function of angle of incidence 6in 
relative to the mean capillary axis. Full symbols: present CTT, open symbols: experimental data: 
7 keV Ne"^"^ (Vikor et al. Experimental transmission rates have been normalized to CTT 

results at 6*^ = 0°. 



B. Angular distributions 

The two-dimensional angular distribution of the guided ions displays anisotropic struc- 
tures (Fig. inj in the O^-Oy plane where 6x is the angle relative to the capillary axis in the 
"scattering plane" , the plane spanned by the incident velocity vector and the capillary axis 
while 6y is the angle perpendicular to the scattering plane. The 6y distribution normal to the 
plane of incidence remains almost constant for all angles. Parallel to the plane of incidence 
{6x) broadening and displacement of the peak from the center of the distribution is found. 
This is in agreement with experiments which show a small deviation of the centroid of the 

n 

scattering distributions towards larger scattering angles for large incidence angles 6in p . At 



ines in Fig. E)) become 
. From the simulation 



certain incidence angles (e.g. 3°) structures in the distribution (soUd 
visible which qualitatively agrees with recent experimental findings 9 
we can delineate the origin of the second maximum as being due to the formation of a small 
secondary charge patch close to the exit surface which deflects projectiles passing by at close 
distance. We point out that double peak structures are stable under variation of Dg and the 
capillary length L although they appear at different angles Om- It remains to be investigated 
to what extent double peaks in the data are non-accidental and can be taken as evidence 
for patch formation near the exit surface. 
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FIG. 6: Two-dimensional angular distributions of 3 keV Ne^^ ions transmitted through PET cap- 
illaries for 6 in = 0° ,1°, 3°, and 5°. The 9^ and 9y directions are defined parallel and perpendicular 
to the plane of incidence given by the ion beam and the capillary axis, respectively. Solid lines are 
projections of the distributions on the ^^.-axis. 

The mean angular spread determined by projecting the data of Fig. IHl onto the 9^ axis 
can be compared with experimental data (Fig. [7j). All distributions in Fig. [7| have been nor- 
malized to the same maximum value. The simulated spread A6out (FWHM) with Fmean = 
is close to the geometric angle A6'o. Including the defocussing effect due to Fmean, AOout in- 
creases by < 1°. Limited by the DC dielectric break-down strength of mylar, the defocussing 
effect of Fmean OH the total spread A9out should not exceed ^ 2°. An additional source of 
spread is the spread in capillary angles. Measurements by Stolterfoht et al. indicate an 
additional spread of up to A6a — 2° due to imperfection in the target preparation can occur 
[2^. Including these contributions, we find A6nut ~ 2.5° (solid line in Fig. [7j), which is still 
somewhat smaller than the value observed in j^. 

C. Transient charging and discharging effects 

Establishing the dynamical equilibrium necessary for guiding takes a finite time interval 
determined by the incoming ion current density and diffusion speed. Studying the transient 
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FIG. 8: Simulated transmission for 3 keV Ne^'^ ions with 9in = 1° averaged over 5 (dashed line) 
and 25 capillaries (full line). The thick line is given by Eq. ESlwith Tg// ~ 1.5 min. 

behavior as the ion beam is switched on and switched off provides additional information 
on the patch formation causing guiding. Fig. |H1 displays the time-dependent transmission 
probability after the beam is switched on at t = assuming that the nanocapillary was 
initially completely discharged. Note that the latter requirement is difficult to assess and to 
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FIG. 9: Simulated transmission for 3 keV Ne ions with — 1° (filled, circles) and — 5° 
(open circles) after the ion beam has been turned off at time Iq = 0. Data has been normalized to 
the transmission at t = to- Lines given by Eg. 1261 

realize in the experiment. This may explain the observation of different transient behaviors. 
The simulated transient discharge behavior for 3 keV Ne'^^ with 6in = 1° and 6^ = 5° is 
shown in Fig. IHl It is probed by first switching off the ion beam at t = to and switching it 
on again after a time interval 6t = t — ta for a sufficiently short time interval in order not 
to disturb the discharge process. Data points are normalized to the transition rate at t = to 
and represent an average transmission over 10 capillaries. 

The charging and discharging behavior can be described in terms of a sirnple analytic 
rate equation model which is a generalization of the model proposed in jl^l . Il4|. To this 
end we approximate the full surface diffusion including the charge absorption at the nearby 
gold layer by an effective surface discharge time which we estimate, in analogy to the bulk 
discharge time (see Eq. [TT|) , by 

Ts = a^Ds . (19) 

The apparent is expected to depend on the angle of incidence 6'j„. On the one hand, the 
extension of the charge patch from the gold layer at the entrance scales as ~ 1/ tan(0i„) ^ 
l/6in, on the other hand, projectiles with a larger energy component normal to the surface 
E± = E sin^ 9in require a larger amount of total charge deposited on the capillary surface 
to be deflected along the capillary axis. The latter requirement is directly reflected in the 



later onset of transmission for 6in = 5° as compared to the switch-on in the case of 6in = 1° 



15 



depicted in Fig.|Hl Taking these effects into consideration we expect for a weak dependence 
on er^- x<l. 

Tlie balance equation for charge qgif) deposited at the surface of the capillary reads 

qsit) = iin- itr - is (20) 

where jj„ is the incident current, jtr- the transmitted current, and js is the current absorbed 
in the capillary surface. The latter has two contributions, current into the bulk and current 
along the surface eventually reaching the gold layer. 
We set 

Js = ( - + - ) Is (21) 

\Ts nj 

Is - lb frtn\ 

qb = -3b = ■ (22) 

n 

Solving Eq. 1201 for switch-on with initial condition qs{t = 0) = 0, we find 

qs{t) = Te// Uin - Jtr) (l " e~'/^^ff) (23) 

where 



1 1 
— + - 

Ts n 

since <C ti,. For the switch-off at t = to we find 



V/ = - + - ^ (24) 



qs{t) = qs{to)-e-'/^^". (25) 

Solving Eq. for any initial condition qb{to) < qs{to) and using the solution for switch 
off of the rapidly varying function qs{t) leads for the total charge-up of the capillary (or, 
equivalently, the capillary transmission) to 

q(t) = q,{t) + q,{t) = qo ■ [C ■ e"*/^^ + {I - C) ■ e-'/^^ff] (26) 

with C = C{qs{to),Tb,Ts) < 1, which is a slowly decaying function (time constant tk) after 
a fast initial decay with time constant Tg// ~ as also found in experiment Q]- A 
comparison with the simulation shows that the switch-on transmission closely follows Eq. 
1221 It should be noted that the delayed onset of transmission cannot be reproduced by a 
simple rate equations. This reflects the fact that the rate equation does not account for the 
threshold behavior for the onset of the trampoline effect. Transmission during switch-off 
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mirrors the total discharge (Eq. I^Uj) . The simulation shows that the time constant is a 
decreasing function of the incidence angle 9in as expected. Fits to the simulation give values 
of Ts ^ 1.5 min and Tf, > 35 min, their order of magnitude consistent with the parameters 
derived from the macroscopic conductivity data Q. The point to be emphasized is that no 
non-linear processes have to be invoked to account for the discharge characteristics. 



IV. CONCLUSIONS 



We have presented a simulation for ion guiding through insulating nanocapillaries within 
the framework of a mean-field classical transport theory. We combine a microscopic trajec- 
tory simulation with macroscopic material properties for bulk an surface conductivities of 
highly insulating materials (PET, "Mylar"). Projectiles hitting the inner wall of the insulat- 
ing material in the early stage of the irradiation deposit their charge on the capillary surface. 
These charges diffuse along the surface and, eventually, into the bulk due to the small but 
finite surface and bulk conductivities of the insulator. Projectiles entering the capillary at 
a later stage are deflected by the Coulomb field of a self-organized charge patch close to the 
entrance of the capillary. Invoking this scenario we are able to reproduce a range of exper- 
imental findings, e.g., ion guiding even for large incidence angles, the temporal decrease of 
transmission during beam-off times, and, in part, a relatively large angular spread of the 
transmitted beam. We have shown, that these results can be interpreted on the basis of a 
linear model including transport of deposited charges along the surface without resorting to 
freely adjustable parameters. Future investigations should address the dependence on ma- 
terial properties. In particular, measurements for other insulating materials (Si02) should 
provide for benchmark data and critical tests for the theoretical framework presented here. 
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